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ABSTRACT 

We present radio observations of GRB 020405 starting 1.2 days after the burst, 
which reveal a rapidly-fading "radio flare". Based on its temporal and spectral 
properties, we interpret the radio flare as emission from the reverse shock. This 
scenario rules out a circumburst medium with a radial density profile p oc r~ 2 
expected around a mass-losing massive star, since in that case the reverse shock 
emission decays on the timescale of the burst duration t ~ 10 2 s. Using published 
optical and X-ray data, along with the radio data presented here, we further show 
that a self-consistent model requires collimated ejecta with an opening angle, 
9j ~ 6° (tj « 0.95 day). As a consequence of the early jet break, the late-time 
(t > 10 day) emission measured with the Hubble Space Telescope significantly 
deviates from an extrapolation of the early, ground-based data. This, along with 
an unusually red spectrum, F u oc z/~ 3,9 , strengthens the case for a supernova that 
exploded at about the same time as GRB 020405, thus pointing to a massive 
stellar progenitor for this burst. This is the first clear association of a massive 
progenitor with a uniform medium, indicating that a p oc r~ 2 profile is not a 
required signature, and in fact may not be present on the lengthscales probed by 
the afterglow in the majority of bursts. 
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Subject headings: gamma-rays:bursts — stars:mass loss 

1. Introduction 

Over the past few years several indirect lines of evidence have emerged in favor of 
massive stars as the progenitors of long-duration gamma-ray bursts (GRBs). These include 
the distribution of offsets of GRBs from their host centers (Bloom, Kulkarni & Djorgovski 
2002), the absence of optical afterglows from dark GRBs (Djorgovski et al. 2001; Piro et 
al. 2002), high column densities toward several GRBs (Galama & Wijers 2001), and the 
inference of very high star formation rates in several GRB host galaxies (Berger, Kulkarni 
& Frail 2001; Frail et al. 2002; Berger et al. 2002). 

Perhaps the most convincing evidence comes from the detection of late-time (t ~ 20 
days) red bumps dominating the optical emission from several afterglows (e.g. Bloom et al. 
1999). While several interpretations of these bumps have been suggested (Esin & Blandford 
2000; Waxman & Draine 2000; Ramirez-Ruiz et al. 2001; Reichart 2001), the preponder- 
ance of spectral and temporal evidence (e.g. Bloom et al. 2002; Garnavich et al. 2002) 
indicates that these bumps are due to radioactive Nickel emission from supernovae (SNe) 
accompanying the bursts. These observations lend support to the collapsar model (Woosley 
1993; MacFadyen & Woosley 1999), in which the core of a massive star collapses to a black 
hole, which then accretes matter and powers the GRB, while the rest of the star produces a 
supernova. 

A seemingly unavoidable consequence of this scenario is that the GRB ejecta should 
expand into a medium modified by mass loss from the progenitor star. To first order, 
the expected density profile is p oc r -2 , arising from constant mass loss rate and wind 
velocity. Extensive efforts have been made to find evidence for such a density profile based 
on broad-band observations of the afterglow emission (e.g. Chevalier & Li 2000; Berger et 
al. 2000; Panaitescu & Kumar 2002). Unfortunately, these studies have been inconclusive 
in distinguishing between a wind density profile and a medium with uniform density, due 
in part to the lack of early observations (particularly in the radio and submillimeter bands) 
and the degeneracy between dust extinction and the intrinsic spectral slope in the optical 
and near-IR bands. Thus, the signature of stellar mass loss remains the elusive missing link 
in the association of GRBs and massive stars. 

To date, the single exception to this disappointing trend is GRB 011121, which provides 
strong evidence for a circumburst medium shaped by a stellar wind (Price et al. 2002), 
and an accompanying SN (Bloom et al. 2002; Garnavich et al. 2002). The reason for 
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these unambiguous results is the combination of exquisite Hubble Space Telescope (HST) 
observations, extensive near-IR data, and dual-band radio data. 

More recently, Price et al. (2002) presented the 7-ray properties and redshift (z = 
0.695 ± 0.005) of GRB 020405, along with multi-band ground-based and HST optical obser- 
vations. The observations between 15 and 65 days after the burst reveal a red bump (with a 
spectrum F v oc z/~ 3 ' 9 ) brighter than an extrapolation of the early data. Price et al. (2002) 
interpret this emission as coming from a SN accompanying the burst, but note that the 
statistical significance of this result depends on the degree of collimation of the GRB ejecta. 
This is because a more collimated outflow results in an earlier steepening of the afterglow 
lightcurves, and hence a more significant deviation at late time. 

In this paper we present radio observations of GRB 020405, which point to a uniform 
density circumburst medium. We also show that the radio, optical, and X-ray data require 
an early jet break, which significantly strengthens the SN interpretation for the late-time 
emission. Combining these two results we conclude that a p oc r~ 2 density profile is not a 
required signature of a massive stellar progenitor. 

2. Radio Observations 

Very Large Array (VLA 1 ) observations were initiated 1.2 days after the burst using the 
standard continuum mode with 2 x 50 MHz contiguous bands. A log of all observations is 
given in Table 1. We used the extra-galactic source 3C286 (J1331+305) for flux calibration, 
while the phase was monitored using J1356— 343. The data were reduced and analyzed using 
the Astronomical Image Processing System (Fomalont 1981). 

3. Reverse Shock Emission in the Radio Band 

In Figure 1 we plot the 8.46 GHz lightcurve of GRB 020405, as well as the radio spectrum 
between 1.43 and 8.46 GHz on day 3.3. The early emission is characterized by two important 
features. First, it is brightest (F u « 0.5 mJy) during the first observation (t m 1.2 days), and 
rapidly fades, F v oc t~ L2±a4 , between 1.2 and 5 days. Second, the spectral index between 
1.43 and 8.46 at t « 3.3 days is (3 ia d —0.3 ± 0.3, and similarly at t « 1.2 days f3 ia d < —0.5 
based on the 8.46 and 22.5 GHz data. 



1 The VLA is operated by the National Radio Astronomy Observatory, a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. 
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The rapid fading and negative spectral slope are atypical for emission from the forward 
shock on the timescale of 1 day. In fact, for typical parameters, the radio band, z/ ra( j, lies well 
below the peak of the synchrotron spectrum at early time, v m m 2QQe l ^ 2 _ 2 e 2 _ l El ! !^t^ d ' i ^ 2 GHz 
(Sari & Esin 2001); here e# = 0.01eB,-2 and e e = 0.1e e; _i are the fractions of shock energy 
carried by the magnetic fields and electrons, respectively, and E = 10 52 i?52 is the afterglow 
kinetic energy. In this regime, the spectrum is F v oc z/ 1 / 3 (or F v oc v 2 if z/ rad < u a , the 
synchrotron self-absorption frequency; Sari, Piran & Narayan 1998). As a result, the early 
radio flux is expected to either increase, as t 1 ^ 2 or t 1 (Sari, Piran & Narayan 1998; Chevalier 
& Li 2000), be flat (Chevalier & Li 2000), or decay as r 1/3 (Sari, Piran k Halpern 1999), 
followed by a steep decline as t^ 1 to t~ 2 when v m crosses the radio band. This behavior 
is observed in the most radio afterglows (e.g. Berger et al. 2000). Since the early radio 
emission from GRB 020405 does not follow this general trend, we conclude that it did not 
arise from the forward shock. 

Instead, we interpret the observed emission as coming from a reverse shock (Meszaros 
& Rees 1997; Sari & Piran 1999) plowing back into the relativisitic ejecta. Similar "radio 
flares" have been observed in several afterglows (Kulkarni et al. 1999; Frail et al. 2000; 
Djorgovski et al. 2001), in particular from GRB 990123 where the flare was shown to be the 
low-energy tail of the V ~ 9 mag optical flash observed at t ~ 50 s (Kulkarni et al. 1999; 
Sari & Piran 1999). In all cases in which a radio flare has been observed, the emission had 
similar properties to that from GRB 020405, namely a bright (F v ~ 0.5 mJy) flux measured 
at early time (t ~ 1 day) followed by a rapid decline. 

Hydrodynamical studies of the reverse shock (e.g. Kobayashi & Sari 2000) have shown 
that the radio emission depends on the properties of the ejecta and circumburst medium. 
In particular, the main implication of the radio flare from GRB 020405 is that it effectively 
rules out a circumburst medium with a density profile p cx r~ 2 (hereafter, Wind). Chevalier 
& Li (2000) have shown that in a Wind environment, with typical afterglow parameters, the 
cooling frequency of the reverse shock, z/ Cj rs ~ 4 x 10 8 t sec Hz, is significantly lower than its 
characteristic frequency, z/ m ,Rs ~ 10 19 ^c Hz. Thus, the emission peaks at ^ c ,rs, with a flux, 
F up ~ 16 Jy, independent of time until the reverse shock crosses the shell at t CI = 5(l + z)A 10 
s (Chevalier & Li 2000); here A = 10A 10 light seconds is the initial width of the shell. 
Following the shell crossing, electrons are no longer accelerated, and since the reverse shock 
is highly radiative (z/ c ,rs ^ m ,Rs) the emission decays exponentially. Thus, strong emission 
from the reverse shock at t ~ 1 — 2 days is not expected in a Wind environment, indicating 
that the early radio emission from GRB 020405 requires a circumburst medium with uniform 
density (hereafter, ISM). 

In addition, based on the flat spectral slope between 1.4 and 8.5 GHz measured at t — 3.3 
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days, we conclude that for both the reverse and forward shocks v a < 1.4 GHz. Otherwise, 
the emission from the reverse shock would be severely attenuated over this frequency range 
by the forward shock, -F„ obs = F Ufim e~ T " , resulting in a significantly steeper spectrum; here 
r„ is the synchrotron optical depth. 

4. Uniform Density Models for the Afterglow Emission 

Using the conditions inferred in §3 we model the radio, optical, and X-ray data with a 
model describing self-consistently the time evolution of the forward and reverse shocks in a 
uniform density medium. We consider the optical data only at t < 10 days since the emission 
at later times is dominated by a much redder component, possibly a SN (Price et al. 2002). 
We return to this point in §5. 

The time evolution of the reverse shock spectrum (F u oc z/ 1 / 3 for v < z/ m ,RS an d F v oc 
u -( P -i)/2 for y > Um R ^ ig described by „ m RS ^ ^(8+5^/7(1+2^ f ^ rs K r (i2+ii 5 )/7(i+2 fl)) 

and the time of peak emission, t p = max[tdur/(l + z),td ec \ (Kobayashi & Sari 2000); here 
tdur = 60 sis the duration of GRB 020405 (Price et al. 2002), t dcc = (3 J B/327rr^n m p c 2 ) 1/3 , T 
is the initial Lorentz factor, and no is the circumburst density. The parameter 3/2 < g < 7/2 
describes the evolution of the reverse shock Lorentz factor, T oc r~ 9 , and the limits correspond 
to adiabatic expansion (g — 3/2) and pressure equilibrium between the forward and reverse 
shocks (g = 7/2). To evaluate t p , z-Vrs (£ p ) , and i^o.Rsfe) we use the physical parameters of 
the ejecta and circumburst medium as inferred from the forward shock emission (see below 
and Table 2), in conjunction with equations 7-9 of Kobayashi (2000) for the thick shell 
case (i.e. when the reverse shock is relativistic and effectively decelerates the shell), and 
equations 15-17 for the thin shell case (i.e. when the reverse shock cannot decelerate the 
shell effectively). We set the nominal values of g = 3/2 and 7/2 for the thin and thick shell 
cases, respectively. Thus, the only free parameter of the reverse shock emission is To- 

For the forward shock we use the time evolution of the synchrotron spectrum in the 
appropriate regime (i.e. spherical ISM for t < tf Sari, Piran & Narayan 1998, and an 
expanding jet for t > tf Sari, Piran & Halpern 1999); here tj, the jet break time, is the 
epoch at which T ~ Qj 1 , and 6j is the half opening angle of the jet. To account for possible 
extinction within the host galaxy, A\? st , we use the parametric extinction curves of Cardelli, 
Clayton & Mathis (1989) and Fitzpatrick & Massa (1988), along with the interpolation 
calculated by Reichart (2001). 

The results of the two models (thin and thick shell cases) are shown in Figures 1 and 2, 
and summarized in Table 2. We find that both models provide an equally adequate fit (with 
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Xmin ~ 1 P er degree of freedom), but r is not well constrained, T ~ 60 — 5 x 10 3 (Figure 3). 

More importantly, both models require collimated ejecta with tj « 0.95 day (hereafter, 
Jet). Models with a significantly wider collimation angle have Xmin ~ 10 P er degree of 
freedom, primarily because they underestimate the flux in the X-ray band by a factor of 
about 20, and cannot explain the radio and optical emission simultaneously. 

At the same time, the Jet models underestimate the radio flux at t > 15 days by about 
4o" on day 16.3, and about 2.5a on day 32. This is due to an apparent brightening of the 

1/2 

radio emission on this timescale. Since in the radio band F v oc n , one possible explanation 
for the brightening is that the forward shock encounters a density enhancement; a density 
increase by a factor of ten is required. The flux in the optical bands would remain largely 
unaffected since v opt < u c (see Table 2), in which case the flux is independent of density. 

Using the parameters of the forward shock emission we calculate E 2 x 10 51 erg, 
e e 0.1, tB ~ 0.5, and n ps 10~ 2 cm" 3 using equations 4.13-4.16 of Sari & Esin (2001). 
The opening angle of the jet is 9j ~ 0.1 (Frail et al. 2001). Using this value we find a 
beaming corrected 7-ray emergy, E^ ~ 3.7 x 10 50 erg (Price et al. 2002), typical for long- 
duration GRBs (Frail et al. 2001), and a beaming corrected kinetic energy, Ek ~ 10 49 erg, 
lower than the typical inferred values of 10 50 to 3 x 10 51 erg (Panaitescu & Kumar 2002). 

5. Implications for the Progenitor of GRB 020405 and GRB Circumburst 

Environments 

In the previous section we did not consider the optical emission at t > 10 days since it 
has a distinct spectrum, F v cx j3~ 3 ' 9±0A , compared to F v oc / g- 1 - 15±0 - 07 for the early afterglow 
data. Moreover, the predicted brightness of the afterglow at late time is lower by a factor of 
about 3 — 5 in the R and I bands compared to the flux measured with HST (Figure 2). These 
two observations indicate that the late-time emission comes from a separate component. 

Price et al. (2002) interpreted this emission as coming from a supernova that occured 
at about the same time as the burst. However, they note that the significance of this 
conclusion depends sensitively on the time of the jet break. Based solely on the optical 
data, these authors were unable to significantly constrain tj. However, our combined radio, 
optical, and X-ray model with tj ps 0.95 day indicates that the SN interpretation is secure. 

We gain further confidence about this interpretation by comparing the late-time emission 
to the optical emission from the Type Ic SN1998bw (Galama et al. 1998). In Figure 2 we 
plot the combined emission from the afterglow of GRB 020405 and SN 1998bw redshifted 
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to z = 0.695. We correct the SN lightcurves for Galactic extinction, E(B — V) — 0.054 
mag (Schlegel, Finkbeiner & Davis 1998), as well as extinction within the host galaxy as 
determined in §4. While SN 1998bw does not provide a perfect fit to the data, the level of 
agreement is remarkable given that there are other effects at play (e.g. an earlier decay for 
a fainter luminosity; Iwamoto et al. 1998). 

The conclusion that GRB 020405 was accompanied by a supernova indicates that the 
progenitor must have been a massive star. However, we also demonstrated that the circum- 
burst density profile is uniform, at least in the range r ~ 10 15 to 10 17 cm. This is the first 
case in which a massive progenitor and a uniform ambient medium have been associated 
directly, leading us to conclude that a strong density gradient, p oc r~ 2 , on the lengthscales 
probed by the afterglow is not a required signature of a massive stellar progenitor. In fact, 
this may explain why the signature of stellar mass loss has not been observed in the majority 
of afterglows. 

The uniform medium around the progenitor of GRB 020405 does not imply that the 
progenitor did not lose mass. In fact, it has been previously suggested that a relatively 
uniform medium can occur downstream from the wind termination shock (Ramirez-Ruiz et 
al. 2001). This region can extend out to 10 18 cm. In addition, a density enhancement is 
expected at the termination shock, which can explain the increased radio flux at t > 15 days 
compared to the Jet model predictions (§4). 

6. Conclusions 

We showed that the early radio emission from GRB 020405 was dominated by the reverse 
shock, and that this directly implies a uniform circumburst medium. The same conclusion 
holds for all bursts in which radio flares have been detected on the timescale of ~ 1 day. The 
broad-band data indicate that the ejecta underwent a jet break at t ~ 0.95 day (i.e. 9j ~ 6°), 
resulting in a significant deviation of the late-time optical emission measured with HST from 
an extrapolation of the model. Combined with the spectral properties of the late emission, 
and a reasonable agreement with the optical emission of SN 1998bw, this indicates that 
GRB 020405 was accompanied by a supernova. Thus, the progenitor of this burst was a 
massive star. 

The association of a massive stellar progenitor with a uniform circumburst medium 
indicates that the tedious search for a p oc r~ 2 density profile may have been partly in vain. 
It appears that a Wind profile does not necessarily accompany every GRB, and in fact may 
not be the case for the majority of bursts. This result, in conjunction with the inferred Wind 
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medium for GRB 011121 which was also accompanied by a supernova, indicates that the 
circumburst environments of GRBs are more diverse than the simple assumption of constant 
mass loss rate and wind velocity; the interaction of the wind with the local environment may 
play a significant role. Deeper insight into the structure of the circumburst medium requires 
rapid localizations and dense, multi-band follow-up. 

EB thanks R. Chevalier and R. Sari for valuable discussions. We acknowledge support 
from NSF and NASA grants. 
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Table I. VLA Radio Observations of GRB 020405 



Epoch 


At 


^0 


Flux Density 


(UT) 


(days) 


(GHz) 


(a* j y) 


Apr 6.22 


1.19 


8.46 


481 ± 36 


Apr 6.25 


1.22 


22.5 


< 300 


Apr 8.36 


3.33 


8.46 


157 ±43 


Apr 8.38 


3.35 


4.86 


160 ±51 


Apr 8.40 


3.37 


1.43 


234 ± 77 


Apr 8.42 


3.39 


22.5 


< 600 


Apr 9.40 


4.37 


8.46 


121 ±39 


Apr 10.28 


5.25 


8.46 


83 ±18 


Apr 13.24 


8.21 


8.46 


47 ±27 


Apr 21.28 


16.3 


8.46 


115 ±21 


May 7.20 


32.2 


8.46 


64 ±20 


May 21.12 


46.1 


8.46 


61 ±26 


May 24.13 


49.1 


8.46 


-9 ±20 


May 27.17 


52.1 


8.46 


-11 ±25 


May 21.12-27.17 


46.1-52.1 


8.46 


0± 14 



Note. — The columns are (left to right), (1) UT date of 
each observation, (2) time since the burst, (3) observing 
frequency, and (4) flux density at the position of the radio 
transient with the rms noise calculated from each image. 
The last row gives the flux density at 8.46 GHz from a 
co-added map of the data from May 21-27. 
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Table 2. Uniform Density Jet Models for GRB 020405 



Parameter 


Thin Shell 


Thick Shell 


Va (HZ) 


10 9 


10 9 


v m [tiz) 


i.y x iu 


j f> v 1 ni2 
4.U X 1U 


v c (Hz) 


3.5 x 10 14 


1.0 x 10 14 


Fujo (A*Jy) 


740 


760 


V 


1.78 


1.73 


tj (days) 


0.99 


0.91 


A^ ost (mag) 


0.28 


0.25 


F a 
1 


1.2 x 10 3 


3.2 x 10 3 


Xmin/ dof 


61.2/61 


58.7/61 


-^iso,52 


0.3 


0.2 


n 


0.05 


0.08 


e e 


0.1 


0.1 




0.3 


0.7 



Note. - Best-fit synchrotron parame- 
ters, and the inferred physical parameters 
for the reverse and forward shock model 
described in §4. The quoted values for the 
forward shock are at t = tj. We give the 
results for the thin shell (with g = 3/2) 
and thick shell (with g = 7/2) cases. In 
both models with fix the value of u a at 1 
GHz based on the flat spectrum in the ra- 
dio band between 1.4 and 8.5 GHz (§3). a 
T is not well constrained and can take a 
wide range of values (Figure 3). 
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Time Since the Burst (days) 



Fig. 1. — Radio lightcurve at 8.46, and the spectrum between 1.4 and 8.5 GHz on day 
3.3 (inset). The solid lines represent the best-fit combined emission from the reverse and 
forward shock model described in §4 for the thin and thick shell cases. The dashed (thin 
shell) and dash-dotted (thick shell) lines show the contributions of the reverse and forward 
shocks separately. 
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Time Since the Burst (days) 



Fig. 2. — Optical and X-ray lightcurves of GRB 020405. The triangles are upper limits, 
while the open circles are late-time observations, dominated by emission from a supernova. 
The solid lines represent the best-fit reverse and forward shock model described in §4, and 
the dashed lines are the combined afterglow flux and flux from SN 1998bw redshifted to 
z = 0.695. The overall agreement indicates that the excess emission is due to a SN similar 
to SN 1998bw. 
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Fig. 3. — Contours (1, 2, and 3 a) of T versus g based on the best-fit models for the reverse 
and forward shock emission in the thin and thick shell cases (Table 2). Both T and g 
are not well constrained. In particular, T has two sets of minima in the thin shell case, 
corresponding to the cases when t p = tdec and t p = idur/(l + z). The range of T values 
in both models encompasses reasonable values based on expectations from pair production 
opacity. 



